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SYNOPSIS 

During the period of June 1982 to March 1983, there was an extensive 
shutdown of the smelters in Sudbury. This report is written to address two ques- 
tions: 

1. Is the impact of the smelter shutdown on wet and dry deposition observable? 

2. Can the smelter impact be quantified? 

Monitoring results from the Acidic Precipitation in Ontario Study 
(APIOS) event and cumulative deposition networks during July to March of 1980-81, 
81-82, and 82-83 were examined. 

Due to the lower time resolution of the cumulative data and the large 
meteorological variability in the data, comparison of the overall averaged results 
of the operating and shutdown periods obtained within a radius of 250 km of the 
smelters failed to yield any evidence of systematic and definitive changes in 
average precipitation concentration due to the smelter shutdown, except for 
copper. Though not quantitative, a comparison of precipitation concentration con- 
tours obtained over the operating and shutdown periods indicated evidence of a 
smelter impact on SO^= and Cu. However, the comparison of the corresponding 
deposition contours did not. 

By classifying the event data by their synoptic meteorological and air 
back trajectory patterns, the potential smelter contribution to wet and dry 
deposition was quantitatively assessed, and is summarized in the following table 
according to site location and parameter of interest. The values shown in this 
table are thought to be near the upper limits of the possible smelter contribution 
(i.e. the actual contribution is less than these values). 



Site 


Parameter 


Sudbury Smelter 


Sudbury Smelter 






Contribution to 


Contribution to 






Wet Deposition 


Dry Deposition 






(% of total) 


(% of total) 


Dorset 


SO<,= 


11 - 12% 


18- 28% 




SO2 


** 


31 - 37% 


Railton 


SOti= 


2% 


♦ 




SO2 


** 


27% 


Longwoods 


SO^= 


* 


* 




SO 2 


»* 


5% 


Fernberg 


scv 


* 


* 




so 2 


** 


4% 


Kapuskasing 


so** 


* 


46% 




so 2 


** 


47% 


Chalk River 


so^= 


12% 


17- 20% 




so 2 


»♦ 


25 - 26% 



* Not estimated due to limited data. 
** Not measured in precipitation samples. 



It should be noted that the stratified analysis was based on only a 
limited number of samples, therefore the estimated percentage contribution to wet 
and dry deposition due to the Sudbury smelters should be considered as tentative 
and should be confirmed when future operating and shutdown data become avail- 
able. 
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1. INTRODUCTION 

Sudbury is an area of major environmental interest because it is the 
location of two smelting/mining companies, namely, INCO Limited and Falcon- 
bridge Limited. Both smelters emit large quantities of SO2, which accounted for 
U5% and 7% respectively of the Ontario SO2 emissions in 1981 (Yap, 1984). 

During the period of June 1982 to March 1983, there was an extensive 
shutdown of the smelters in Sudbury. Even though during parts of this period, there 
were SO2 emissions associated with the Falconbridge operation (June 1 to 26, 1982 
and January to March 1983), these emissions were considered to be small relative 
to those of the regular total Sudbury Smelter emissions. In order to relate receptor 
observations to source emissions in the acid rain debate, it is pertinent to examine 
the empirical data to address the following two questions: 

1 . Is there any observable effect on wet and dry acid deposition as a result 
of the Sudbury smelter shutdown? 

2. Can the extent of impact be quantified geographically? 

This report is a summary of some attempts made to answer these ques- 
tions using the wet and dry deposition data from the Acidic Precipitation in On- 
tario Study (APIOS) networks, Environment Canada's Air and Precipitation network 
(APN) and Ontario Hydro network. 
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2. DEPOSITION NETWORKS AND DATA 

The Air Resources Branch of the Ontario Ministry of the Environment 
has been quite active in the monitoring of precipitation quality in the Sudbury area 
since the late 70*s (Chan, 1982). It was recognized then that acid rain was not a 
problem of a local origin, but rather it was related to long range transport of air 
pollutants. Hence province-wide networks were established in 1980 to monitor 
both wet and dry deposition. 

There are two networks operated by the Air Resources Branch in the 
province of Ontario: a cumulative network which samples on a 28 day basis (this 
network sampled on a monthly basis from network inception till January 5, 1982) 
and an event network which samples on a daily basis. In both networks, wet 
deposition is directly measured and dry deposition is inferred from ambient air 
concentration measurements. The primary function of the cumulative and event 
networks is respectively to quantify the long term deposition pattern and to 
determine the origin of the pollutants being measured. At this moment, the 
cumulative network has 36 wet monitoring stations of which 23 are equipped with 
air measurement instrumentation. The event network consists of 15 wet deposition 
sites and 4 of them also sample air on a daily basis. The spatial distribution of the 
stations of which results are being used in this analysis is shown in Figures 2-1 and 
2-2 and details of the networks are documented elsewhere (Chan et. al, 1982 and 
1984). 

Chemical analyses of the collected samples are carried out in the 
Laboratory Services and Applied Research Branch (LSARB). Laboratory and field 
data are merged and stored in the Ministry's Sample Information System (SIS). The 
laboratory data are further subject to stringent validation/screening criteria such 
as sample integrity checks (e.g. ionic balance, theoretical vs observed 
pH/conductance, free hydrogen vs total hydrogen ions etc), gross limit checks and 
Dixon Ratio test (Kirk 1983 and 1984). Outliers are flagged but not deleted and 
are retained in the data base together with the validated data. 
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3. DATA ANALYSIS 

The analysis performed here places a higher emphasis on the wet 
deposition results because there are more wet deposition data available than dry 
deposition/air concentration data. 

In order to fully understand the relationship between source and 
receptor, information on source emissions, intermediate processes such as 
dispersion/transport, chemical transformation, and deposition at the receptor must 
be available. One difficulty often encountered in source impact analysis is the 
large variation in the data due to meteorology. This point can be illustrated by the 
outcome of the gross averaging of the observed results for periods corresponding to 
various seasons and source emission configurations, and is the theme in Section 4. 

An examination of the gradient of concentration/deposition as a fun- 
ction of distance from the source yields qualitative information on the potential 
impact of the Sudbury sources. Especially by taking ratios of data to a tracer such 
as NC>3 _ which is not emitted by the smelters, the variability due to meteorology is 
partially accounted for. This approach is described in Section 5, using the 
cumulative wet deposition data. 

One approach of examining the source/receptor relationship is by 
mathematical modeling, holding the meteorological component constant but adjust- 
ing the emissions factor (Ellenton and Misra 1984). An alternative approach is to 
group the results under similar meteorology and compare different emission 
configurations. A first-cut attempt at this is by simple air parcel trajectory 
analysis (Kurtz and Yap, 1984). 

Due to the uncertainties associated with simple air parcel trajectory 
groupings, data need to be examined further with a view to the concurrent synoptic 
pattern, in order to estimate quantitatively the extent of the impact by the point 
source (in this case Sudbury) on a receptor (in this case monitoring site such as 
Dorset). Sections 6 (5c 7 summarize results using this detailed stratification scheme 
to estimate the impact on wet and dry deposition of the Sudbury source on distant 
sites. 
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4. STATISTICAL ANALYSIS TO EVALUATE SUDBURY SHUTDOWN IMPACT 
USING CUMULATIVE NETWORK DATA 

4.1 Analysis: 

Results obtained at stations from the APIOS cumulative network within 
a radius of 250 km from Sudbury are grouped together on a seasonal basis (e.g. 
January to March, April to June, July to September and October to December and 
these periods are labelled in the figures as 02, 05, 08, 1 1 respectively). 

Table 4-1 summarizes the observation periods, number of sample points 
used, precipitation-weighted mean concentrations (average concentration weighted 
by the precipitation depth) and the number of missing data points in this 
calculation. Of particular interest are the 3 shutdown periods (8208, 8211 and 
8302) and their comparisons with respect to the preceeding operating periods. The 
results for selected parameters are represented in Figures 4-1 to 4-5 for easy 
visual comparison. The reason for grouping data in the reported fashion is to allow 
a comparison for the same seasons, and thus decrease variability due to season- 
ality. Table 4-2 shows whether or not there was a statistically significant decrease 
in concentrations during the shutdown period. 

4. 1 . 1 . Hf± (free hydrogen ion); Figure 4-1 

Seasonality is noted in the data; highest values occur in the summer or 
in the fall. The effect associated with the smelters' shutdown on Hf + is somewhat 
irregular; Hf+ decreased in the 08 period but increased in the 11 and 02 periods. 
The observation may reflect an irregular neutralization pattern due to the presence 
of alkaline materials (e.g. NH. , Ca , etc.), suggesting that H f + is perhaps not a 
good parameter to be used for source impact evaluation. 
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4.1.2.H * (total hydrogen ion); F igure fr-2 

The trend in H is again quite irregular for the same reason given 
above for H* + and therefore should not be used as a definitive indicator of smelter 
impact. 

fr.l.3.SOft=; F igure 4-3 

Except in the case of the 08 period where a consistent decrease was 
observed during the shutdown with respect to the operating periods, both the 1 1 
and 02 periods exhibit a high-low-high pattern which is not easily interpretable. 

4.1.4. N-NQ 3 "; F igure 4-4 

Seasonality is again observed of this parameter. Even though there is a 
mild variation for each of the seasons (08, 11, 02) between the shutdown and 
operating periods, the change is not statistically significant, being consistent with 
the fact that only negligible emissions of NO originate from the Sudbury smelters. 

4.1.5.Cu: F igure 4-5 

Cu is a parameter characteristic of smelter operation. In the display of 
the average seasonal concentration, there is a consistent decrease showing the 
effect of smelter shutdown on copper concentrations even though the decrease is 
not always statistically significant. 

4.2. Summary 

An attempt was made to examine whether the effect, due to the 
shutdown of the Sudbury smelters on precipitation quality, is observable. The 
grouping of data from stations within a 250 km radius of Sudbury suggests no 
systematic, statistically significant decrease during the shutdown period, except 
for Cu. This points to the fact that because of the large variability in the results 
associated with meteorology and other factors, more sophisticated approaches such 
as those using stratification schemes have to be used. 
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5. IMPACT OF SUDBURY EMISSIONS ON WET DEPOSITION IN ONTARIO 

5.1 Approach 

This analysis attempt employs wet deposition data from the APIOS 
cumulative sampling network. Precipitation depth weighted concentration 
and wet deposition contours of selected parameters were calculated for both 
the shutdown period and the corresponding operating periods in the 
preceeding years. The period selected is from July to March because of data 
availability. Only stations with at least five (for 80/81 and 81/82 periods) and 
six (for 82/83 period) data points respectively were retained in the contour 
calculations. The calculation of the weighted concentration is straight- 
forward, by using the available data on precipitation gauge depths and 
observed concentrations. Deposition was calculated by multiplying the 
weighted concentration with the corresponding total precipitation depth for 
the period of interest. Because of the complexity in the interpretation of 
deposition data (which depends not only on concentration, but also on 
precipitation amount, which can be quite variable), more emphasis will be 
placed on the discussion of the concentration contours. 

5.2. Concentration Contours 

Plots for Hf + (free hydrogen ion), SO$ s , N-NO3- and Cu are discussed 
below. Hue to the fact that the free hydrogen ion concentration is influenced 
by potential neutralization by alkaline materials (e.g. NHz, + , Ca ++ , Mg++) 
which may be both of local or long distance origins, these plots do not display 
readily recognizable patterns from which definitive conclusions regarding 
emissions impact can be drawn. 

The following summarizes patterns for the three July to March periods 
of 1980/81, 81/82, and 82/83 of SO^=, NO3-, and Cu. Missing data are 
replaced by "-9.00" in the figures. Station identification is summarized in 
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Table 5-1. Aside from the regular APIOS data, results obtained at stations 
near Sudbury (//37 Burwash, #38 Lively and //39 Hammer) which are retained 
from the Sudbury Environmental Study are also used in the contours 
production. 

5.2.1.SO(£ 

Figures 5-1 to 5-3 are precipitation depth weighted concentration 
contours for the three periods of interest. Common features are: 
(i) there is a general S to N concentration gradient, decreasing from S 

toN. 
(ii) there is an elevation around the vicinity of Sudbury - the elevation 
during the shutdown period seems to be less pronounced with respect to 
that of the operating periods and may be attributed to possible effects 
of local dust from the tailings areas near the smelters. 

Both these facts together qualitatively point to the effect of the 
Sudbury emissions on SO^= precipitation concentration. To show the relative 
changes during the shutdown period with respect to the operating periods, 
ratio plots (Figures 5-4 and 5-5) are generated. Except for sites which are 
located at the extreme north and south of the province, ratios of the 
concentrations between the shutdown period and the operating periods are all 
less than unity. Some stations (//23 in Figure 5-4 and //23, 37, 21, and 12 in 
Figure 5-5) along the east side of Georgian Bay have anomalously high ratios 
and the reason for this is unclear. 

5.2.2.N-NO31 

Concentration plots for N-NO3- are shown in Figures 5-6 to 5-8 for the 
three periods. Again a negative S to N gradient is observed, but the elevated 
concentration near Sudbury as observed in the SO^= plots is not seen here. 
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This is consistent with the fact that little NO x is emitted from the Sudbury 
sources. The elevated concentration at station 23 (Killarney) is persistent for 
all three periods but is not easily interpretable. 

Figures 5-9 to 5-10 represent ratio plots of N-NO3- concentrations for 
the shutdown vs. the operating periods. It is of interest to note that, except 
in the vicinity of Sudbury, ratios are in general less than unity on the order of 
0.7 to 0.9. The elevation in the Sudbury area is not easily explained, but 
perhaps can be related to meteorological conditions which are specific to the 
area. 

5.2.3. SO ft Z/N-NO ^: 

In view of the fact that both NO3" and SO^= decrease in concentration 
to a certain extent during the shutdown period and that the decrease may be 
attributed to both emissions and meteorological effects, it is desirable to 
examine ratio plots of SO^= to N-NO3" (Figures 5-11 to 5-13). One important 
feature of this approach is that it partially eliminates the effect of 
meteorology because both SO^ = and N-N03" experience similar meteor- 
ological conditions during the same period. In that sense, N-NO3- (which is 
not emitted from the Sudbury smelters) can be used as a tracer to examine 
the impact of Sudbury Smelter emissions on SO^ = concentration. The 
relatively higher SO^~/N-NOy ratios around Sudbury with respect to other 
parts of Ontario during 80/81 and 81/82 but not in 82/83 indicate local 
impact. 

By taking ratios of the ratios of SOi^ = /N-NOy for the shutdown period 
and operating periods, Sudbury emission impact on precipitation quality 
becomes readily discernible. In Figures 5-14 to 5-15, except in the southern 
and northern most parts of Ontario, these ratios are less than unity with a 
minimum (about 0.7) centered around the Sudbury area. The "unity" contour 
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line should be intrepreted as a boundary beyond which impact due to Sudbury 
sources in not readily discernible because of variability/uncertainty 
associated with meteorology and sampling. A further point to note is that 
the unity contour differs in the two comparison years, suggesting that the 
results of the Sudbury smelter impact should be regarded as qualitative only. 

5.2.4.CU 

Cu, being one of the major emittants characteristic of Cu-Ni smelter 
operation, is expected to show some definite pattern manifesting smelter 
impact on precipitation quality. The approach taken here is similar to that 
taken for SO^=, i.e. plotting out Cu concentrations for the three periods of 
interest and ratioing Cu to a tracer (N-NO3"). 

Figures 5-16 to 5-18 show that a maximum exists around Sudbury with 
general decreases in locations away from the smelting centre. The mild 
maximum during the shutdown period is possibly a result of tailings in the 
surroundings of the smelters. Taking ratios of results of the shutdown period 
to those of the operating periods leads to contours of less than unity ratio 
with a minimum around Sudbury, suggesting smelter effects. Plots of Cu-to- 
N03~ ratios which are not included here exhibit a similar pattern to that of 
the SO^-to-N03"-plots, again pointing to the observed effects of Sudbury 
smelters. 

However, it is somewhat puzzling that ratios in these plots are almost 
all less than unity across the province and that Ni concentration patterns 
similar to those of Cu are not observed. Considering that Cu and Ni are both 
major emittants from the smelters and that they both fall in a similar size 
range (i.e. similar transport and scavenging patterns are to be expected), the 
fact that the pattern exists only for Cu and not for Ni, and also the general 
decrease in the Cu-to-N03 - ratio across the province, are observations 
difficult to reconcile. 
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5.3. Deposition Contours 

Deposition contours of sulfate are given in Figures 5-19 to 5-21 for the 
three periods and those of their ratios are given in Figrues 5-22 to 5-23. 
Though the pattern is not as regular as that in the concentration case, 
nevertheless a decrease in deposition is still observed in the vicinity of 
Sudbury showing the impact of smelter emissions on wet deposition. Results 
at station #21 in 1980/81 and at station //23 in 1981/82 and 1982/83 may be 
anomalous. Similar deposition and deposition ratio plots for N-NO3- are 
shown in Figures 5-24 to 5-26, and 5-27 to 5-28 respectively. There is in 
general a south to north negative gradient except for stations along the coast 
of Georgian Bay and part of Lake Huron in the periods of 80/81 and 82/83. 
The reason for this is unclear but may be related to meteorological factors. 
Because of this, stations in the neighborhood of Sudbury display larger than 
unity ratios in the deposition ratio plots between the shutdown year and the 
preceeding operating years. 

5.4. Concluding Remarks 

The approach taken here yields qualitative but not quantitative results 
showing Sudbury emissions impact. This is mainly due to the large variability 
in the data associated with factors such as meteorology and sampling. 
Quantitative evaluation of the Sudbury smelter impact on wet deposition can 
only be obtained by detailed stratification of the event network data, which 
is the theme of Sections 6 and 7. 
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6. METEOROLOGICAL ANALYSIS OF DAILY PRECIPITATION DATA 

6.1 Dorset Daily Precipitation Data Analysis 

Simple trajectory analysis has been carried out using several years of 
Dorset precipitation data (Kurtz and Yap, 198*0. In this section, the contribution 
of the Sudbury sources is assessed by carrying out a more detailed meteorological 
stratification. The purpose of this stratification is to minimize the variability 
associated with other emission sources and meteorological factors, and to isolate 
the effect of the Sudbury sources (see Appendix 1 for a discussion of the data 
stratification techniques). 

6.1.1 Data Stratification 

The Sudbury sources have the highest probability to impact the Dorset 
site when the synoptic pattern is one of three types: a post-cold front, pre-ridge 
situation; a cold front passing through the area; or a low pressure system to the 
northeast (Figure 6-1). Since a cold front is formed due to a cold air mass (flowing 
toward southeast) pushing and lifting a warm air mass (flowing toward north or 
northeast), the frontal precipitation samples are affected by sources from south or 
southeast even when the surface geostrophic wind back trajectories are in the west 
or northwest sectors. Although the third synoptic type above does not have this 
problem, the number of available samples is very small. Hence, only the post-cold 
front, pre-high pressure ridge type is selected here. However, this synoptic pattern 
may be associated with many different trajectory patterns. In other words, the 
samples may still be affected by sources from many directions, even for this 
particular synoptic situation. In order to improve the chance of identifying the 
Sudbury smelters 1 effects, we use also trajectory information to stratify the data. 
By using both the synoptic and trajectory patterns, we are able to minimize both 
source variability and meteorological variability, and hence, the between-year 
variation of the stratified data. 
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6.1.2 Criteria of Selection 

According to the above discussion, the data selection criteria are: 

(1) post-cold front, pre-high pressure ridge type synoptic pattern; 

(2) the envelope of 48 hour back-trajectories calculated every 6 hours 
during the 24 hour sampling period should include Sudbury; 

(3) the 48-hour back-trajectories should at all times be within the 90° 
sector between west and north (to minimize the possibility of including 
cases where air parcels with recent history from the high emission 
density areas to the south have entered the northwest sector). 

6.1.3 Data Analysis 

6.1.3.1 Selection of the Statistical Test 

The significance of the differences of the stratified samples and their 
concentration means during the shutdown and the operating periods should be 
examined before using them to estimate the Sudbury smelter shutdown effect. In 
general, except for testing whether or not there is a difference between 
populations, all statistical tests for the means are based on the assumption that the 
samples are random, and most are based on the additional assumption that the 
populations have normal distributions. Due to the fact that the populations of the 
present data do not have normal distributions, and each parameter deviates from 
normality differently, the use of these statistical tests based on a normal 
distribution will lead to wrong conclusions. 

The statistical test of the significance of the means is not necessary if 
two stratified samples have significantly different population distributions and one 
sample has a much wider data range than the other, since this condition will pro- 
vide significantly different means. Since the main purpose of the detailed data 
stratification is to minimize the meteorological variability and to maximize the 
Sudbury smelter effect, then for the present analysis to be successful, the 
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stratified data sets should have distinct data ranges, e.g. the SO^ = concentration of 
the shutdown period should be much lower than that of the operating period. In 

other words, for the successfully stratified data sets, a statistical test of the 
population differences is sufficient for examining the significance of the Sudbury 
smelter contribution. For the stratified data sets which do not have distinct data 
ranges, the evaluation of the Sudbury smelter shutdown effect cannot be per- 
formed. 

There are many non-parametric tests available in the literature for 
testing population differences (see Polland, 1977). Since non-parametric tests are 
robust (not sensitive to slight deviations from the assumed population distribution), 
all of them should provide similar test results. To justify this point, these methods 
have been tested by using the SO^= and NO3- concentrations of the stratified 
precipitation samples collected in the three periods of interest (July 80 -March 81, 
July 81 - March 82, and July 82 - March 83) at Dorset. These methods include the 
Mann-Whitney, Kruskal-Wallis, Median 2-sample, Median one way, Vander Warden, 
A NOVA and Wilcoxon 2-sample tests. All of these tests provide the same results 
namely that: (1) the SO^= concentrations were significantly lower when the 
Sudbury smelters were shut down, but the NO3" concentrations were not 
significantly different between shutdown and operating periods; (2) both SO/p and 
NO3" are not significantly different between the two operating periods. This result 
indicates that any of these methods may be used for testing the population 
differences. In this study, we select the median 2-sample test since it is available 
in the standard 'SAS' statistical package (SAS, 1979). 

6.1.3.2 Significance of the Sudbury Smelter Contribution 

According to the criteria of selection described, the number of 
available cases for the 3 periods of interest (July-March of 80/81, 81/82 and 82/83) 
is 9, 3 and 5 respectively. (During the second period there were actually 5 
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meteorologically acceptable occurrences, but on 2 of these occurrences, precipi- 
tation chemistry data were missing.) The SO^= and N-NO3- concentrations in 
these samples are plotted in Figure 6-2 (a and b). The significance of the 
differences among the samples was tested by using the Median 2-sampIe non- 
parametric test. The results (see Table 6-1) indicate that the SOfy- concentrations 
of the shutdown period are significantly lower than those of the operating periods, 
while the NO3- concentrations of all three periods are not significantly different. 
Since the Sudbury smelter contribution to NO3- is negligible, this observation gives 
credence to the meteorological stratification scheme. The precipitation depth 
weighted mean concentrations of these samples are listed in Table 6-2. The 
precipitation depth-weighted mean SO^= concentrations, for the Sudbury smelter 
operating periods, are comparable, but that for the shutdown period is only about 
31 to 33% of those in the operating periods, indicating that the reduction of SCty = 
concentration in the shutdown period is related to the smelter SO2 emissions. 

The Sudbury smelter shutdown effect cannot be shown if one tries to 
compare directly the deposition data rather than the concentration data. This is 
due to the variability of the precipitation amount, which introduces a large amount 
of additional noise into the data analysis. For purposes of illustration, a 
comparison of deposition similar to that for concentration shown above has been 
made. The results (Figure 6-3 and Table 6-4) indicate that there is no statistically 
significant difference of sulfate and nitrate deposition between the smelter shut- 
down and operating periods (as expected due to the additional "noise" in the 
depositions data). 

6.1.3.3 Estimation of the Sudbury Smelter Contribution to the Wet Deposition 
at Dorset 
The Sudbury contribution to the wet deposition at Dorset for each period of 
interest was estimated by multiplying the difference of the precipitation depth- 
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weighted concentrations during the shutdown and operating periods (Table 6-2) with 
the total precipitation depth accumulated over sampling days during which there 
were trajectories coming from the Sudbury region (that is, only criterion 2 of 
Section 6.1.2 had to be satisfied). A number of points should be noted in 
connection with this calculation: 

1. Since many more cases are included in calculating the total precipitation 
depth associated with the NW quadrant than the limited number of events 
identified by the rigorous stratification scheme of Section 6.1.2, the total 
depth used here is considerably greater than the accumulated values of the 
limited events of Section 6.1.2. 

2. The estimate of Sudbury contribution resulting from this calculation is 
expected to be on the high side, since not every one of the events having a 
surface geostrophic trajectory in the NW qradrant will be associated with 
such a high Sudbury impact as in the vigorously stratified data set. 

3. The reason that the rigid criteria of Section 6.1.2 were relaxed for this 
estimate of Sudbury impact on deposition, is the uncertainty about the proper 
height at which trajectories should be calculated for an elevated source such 
as the INCO smelter. When the plume stays aloft during the precipitation 
event, the 850 mb back-trajectory may be more relevant than the surface 
geostrophic trajectory. In such cases, the 850 mb back-trajectory could be 
over Sudbury, while the surface geostrophic trajectory is from a more 
northerly direction. Hence the sector was "widened" to include all the events 
from the NW quadrant. 

Results of this calculation indicate that the Sudbury smelter contribution to 
the Dorset wet deposition (see Table 6-3) was about 12 and 11% as a probable 
upper-level estimate for the first and second operating periods, respectively. 
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6.2 Railton Daily Precipitation Data Analysis 

6.2.1 Data Stratification 

Due to the fact that the Sudbury smelter contribution to a receptor is 
inversely proportional to the distance between the receptor and the smelters, it is 
expected that the Sudbury smelter contribution to Railton is smaller than that at 
the Dorset site. The chances of identifying the Sudbury smelter effects by 
meteorological stratification are also smaller. This is due to the fact that there 
are more important meteorological factors which have to be included in the 
criteria of data stratification. For study periods of less than a year, as is the case 
here (since there was only one shutdown period) the chances of success of data 
stratification are further limited by data availability. 

6.2.2 Criteria of Selection 

The criteria of data selection for Railton are: 

(1) Post-cold front, pre-high pressure ridge type synoptic pattern (See Figure 6- 
*>, 

(2) The envelope of 48-hour back trajectories calculated every six hours during 
the 24-hour sampling period should include Sudbury; 

(3) The 48-hour back trajectories should at all times be within the 90° sector 
between west and north. 

The above criteria were used to estimate the contribution of Sudbury 
smelter to the precipitation concentrations at Railton by comparing values 
observed during the operating and shutdown periods for the set of samples 
satisfying these criteria. 

As with the Dorset data, an estimate of the contribution of Sudbury to 
the wet deposition at Railton (probably on the high side) was made by using accum- 
ulated depths from all events where the 48-hour back trajectories were, at least 
one time during the sampling period, within the 90° sector between west and north. 
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6.2.3 Data Analysis 

6.2.3.1 Significance of Sudbury Smelter Contribution 

According to the criteria of selection, the number of cases for the 
three periods of interest (July-March of 80-81, 81-82, and 82-83) are 6, 6, and 3, 
respectively. It was found that both the SCty= and N-NO3- concentrations are 
lower in the shutdown period. Since the Sudbury smelters are not a significant 
NO3" source, the decrease in NOy concentration indicates that factors other than 
the smelter shutdown are also responsible for the reduction of concentration during 
the shutdown period. This also suggests that further data stratification is 
necessary. Unfortunately, this was not possible because the sample size is too 
small. 

The significance of these differences is tested by using the median 2- 
sample test. Results are listed in Table 6-1. These results suggest that the 
Sudbury smelter effects are significant for the first operating period, while factors 
other than emissions from smelters are the main contributors for the second 
operating period. In other words, the second operating period data should not be 
used in estimating the Sudbury smelter contribution. 

6.2.3.2 Estimation of Sudbury Smelter Contribution to the Wet Deposition at 
Railton 

The Sudbury smelter contribution to the wet deposition at Railton for 
each period of interest was estimated by using the approach similar to that at 
Dorset (see section 6.1.3.3). The precipitation depth weighted mean concentration 
for each period is given in Table 6-2. These results indicate that the precipitation 
sulfate concentration associated with the most likely Sudbury smelter affected 
samples is about 12% (from a comparison of the 80-81 and 82-83 data), and the 
Sudbury smelter contribution is approximately 2% of the total wet deposition for 
the first operating period (July 1980 - March 1981) (see Table 6-3). 
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6.3 Longwoods Daily Precipitation Data Analysis 

6.3.1 Data Stratification 

The synoptic pattern for the cases most likely to be affected by 
Sudbury is the post cold front pre-high pressure ridge condition when Longwoods is 
located to the west or northwest of a low pressure system (see Figure 6-5). The 
48-hour back trajectories for these cases should be in the north sector (90°). 

6.3.2 Criteria of Selection 

The criteria of selection for the samples at Longwoods are: 

1. post cold front pre-high pressure ridge type synoptic pattern (i.e Longwoods 
is located to the west or northwest of a low pressure system); 

2. the 48-hour back trajectories are in the 90° sector between northwest and 
northeast. 

6.3.3 Data Analysis 

6.3.3.1 Significance of Sudbury Smelter Contribution 

According to the criteria of selection, the number of cases for the 
three period of interest CJuly to March 80-81, 81-82, and 82-83) are 3, 6, and 1, 
respectively. Since there is only one case in the third period, the significance test 
was not carried out for the third period. The results (see Table 6-1) of the median 
2-sample test indicate that there is no significant difference between the first and 
the second periods (the operating periods). The estimation of Sudbury smelter 
contribution to deposition at this site would not be meaningful due to the fact that 
there is only one case during the shutdown period that falls within the criteria of 
selection. 
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6.4 Fernberg Daily Precipitation Data Analysis 

Since the Fernberg daily precipitation sampling began only in October 
1981, the sample size is not large enough for carrying out the detailed meteoro- 
logical stratification. Only synoptic patterns were used for data stratification. 

6.4.1 Data Stratification and Results 

The synoptic patterns corresponding to the most possible Sudbury 
smelter affected cases are (see Figure 6-6): 

1. Fernberg to the northwest of a high pressure system; 

2. Fernberg to the southwest of a high pressure system; 

3. Fernberg to the northeast of a low pressure system. 

There were three cases for the shutdown period, yet only one case for 
the second operating period and zero for the first period, since no samples were 
collected. These data do not indicate any shutdown effect since both SO^ = and N- 
NO3* concentrations are lower in the shutdown period. 

6.5 Kapuskasing Daily Precipitation Data Analysis 

An attempt was made to examine the smelter impact on precipitation 
samples collected from the Ontario Hydro's Kapuskasing site. This site is attrac- 
tive because of its location being north of Sudbury, which is more frequently 
impacted by air masses passing through Sudbury. 

Due to the fact that many precipitation samples of the operating period 
are missing (from October to December as well as some events in the other 
months), detailed stratification of the samples was not possible (Jarv, 1982 and 
Handy, 1984). In total there were only 7 and 16 possible Sudbury smelter-affected 
cases in the shutdown and operating periods, respectively. However, no similar 
patterns could be found for these periods. Hence, the analysis was not performed. 
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6.6 Chalk River Daily Precipitation Data Analysis 

6.6.1 Data Stratification 

Precipitation data (Barrie et al., 1982) obtained from the Evnironment 
Canada APN station site at Chalk River, east of Sudbury, were also assessed for 
Sudbury smelter impact. The synoptic pattern for the most-possible Sudbury 
smelter-affected cases is given in Figure 6-7. Precipitation samples collected in 
the shaded area in Figure 6-7 are the ones most probably affected by the Sudbury 
smelter emissions, since the back-trajectories at Chalk River are then from the 
Sudbury region. 

6.6.2 Criteria of Selection 

The criteria of selection of the for Chalk River are: 

1. north or northwest of a high pressure system (which is located to the south of 
a low). 

2. the 48-hour back trajectories should be at all times within the west sector 
(45°). 

These cases were used to estimate the excess Sudbury smelter 
contribution to the precipitation concentration. The possible Sudbury smelter 
contribution to the wet deposition at Chalk River was estimated by the product of 
the excess concentration and the total precipitation depth, using all cases where 
the 48-hour surface geostrophic wind back trajectories during the sampling period 
were least one time within the 90° sector between northwest and southwest. 

6.6.3 Data Analysis 

6.6.3.1 Significance of the Sudbury Smelter Contribution 

In general, the SO^ = concentrations of the shutdown period are lower 
than those of the second operating period. However, the NO3* concentrations for 
both periods are in the same range. The results of the median 2-sample test (see 
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Table 6-1) indicated that the SO^= concentrations of the shutdown period are 
significantly different from those of the second operating period. However, the 
NO3" concentration of the shutdown period are not significantly different from 
those of the second operating period. There is no significant difference between 
other periods. These results indicate that the Sudbury smelter contribution to the 
sulfate precipitation concentration in the second operating period is significant. 

6.6.3.2 Estimation of the Sudbury Smelter Contribution to the Wet Deposition 

at Chalk River 

The mean SO^= concentration and the percentage difference between 
the shutdown and operating periods are given in Table 6-3. The Sudbury smelter 
contribution to the SO^= concentration for the cases stratified according to the 
above criteria is about 49%. There are 49 possible Sudbury smelter-affected cases 
in the second operating period based on meteorological analysis, however, only 1 1 
samples are available. The total precipitation depth associated with possible 
effects due to Sudbury is estimated by multiplying the average precipitation depth 
of the available 1 1 samples with the total number (49) of possible Sudbury smelter- 
affected cases. The Sudbury smelter contribution to the wet deposition is then 
estimated by multiplying this estimated total possible precipitation depth with the 
excess concentration. Results are given in Table 6-3, indicating that the Sudbury 
smelter contribution to SO^= wet deposition at this site is about 12%. 
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7. METEOROLOGICAL ANALYSIS OF AIR CONCENTRATION DATA 

7.1. Dorset Air Sample Analysis 

7.1.1 Data Stratification 

Due to the fact that the 6 or 1 2 hour back-trajectory sector may not be 
a sufficient indicator of the contributing source and meteorological conditions (see 
Appendix 1), in this section, the contribution to dry deposition from the Sudbury 
sources was examined by employing a more detailed data stratification. As 
discussed in Appendix 1, a specific (48-hour backward) trajectory pattern generally 
is associated with a specific synoptic pattern variation. Thus the trajectory 
pattern was here used to take into account both the contributing sources and 
synoptic pattern variations. 

The trajectory pattern for which the Sudbury sources are expected to 
have the greatest impact at Dorset is the one where trajectories during the 
sampling period are all bounded by the west and northwest directions (see examples 
in Figure 7-1). The second most possible trajectory pattern is that which has the 
48-hour back-trajectories in the NW sector most of the time during the sampling 
day (Figure 7-2). The air concentrations are given in Figures 7-3 to 7-5 for events 
associated with the patterns of this type. From these figures, we can see that 
during the shutdown period the air SO^ = and SO2 concentration were much lower 
than those of the operating periods. However, the NO3" concentrations during the 
shutdown period were higher than those of the second operating period and slightly 
lower than those of the first operating period. The median 2-sample test of these 
samples shows that the shutdown period SO2 and SO^= concentrations are 
significantly different from those of the operating periods (Table 7-1). However, 
the NO3" concentration of the shutdown period are not significantly different from 
the second operating period and significantly different from the first operating 
period. These results indicate that the SO^ = and SO2 concentration difference 
between the shutdown and operating periods may be attributed to a Sudbury 
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smelter contribution during the second period; however, during the first period 
differences may have been caused by both Sudbury and other factors. In other 

words, the Sudbury smelter contribution may be estimated for the second period 

only. Hence, only the analysis for the second period is given in the next section. 

7.1.2. Estimation of the Sudbury Contribution to the Dry Deposition at 

Dorset. 

The contribution of the Sudbury smelters to the dry deposition of SO2 
and SOty = at Dorset was estimated by a method similar to that used for wet 
deposition. Since dry deposition may be estimated by multiplying the appropriate 
dry deposition velocity with the air concentration, and the dry deposition velocity 
may be assumed to be constant, on the average, from year to year for the same 
season, the percentage contribution to dry deposition may be estimated by the 
percentage contribution to the observed concentration (Table 7-3). By comparing 
the shutdown (column 7) and operating (column 5) period concentrations for the 
subset of samples on which the Sudbury sources are expected to have the greatest 
impact, the excess concentration due to the Sudbury source was obtained in column 
8. This concentration was then multiplied with the total number of cases where 
Sudbury could have had a possible impact on Dorset (column 10) and divided by the 
average concentration measured during the operating period times the total 
number of measurements during the operating period. The criterion of selecting all 
the possible Sudbury affected cases is that the 48-hour back trajectories should at 
least one time during the sampling period be within the 90° sector between west 
and north. This criterion takes into account upper level trajectories including all 
the possible Sudbury effects on air concentration through direct transport, and 
indirect mixing from the surface to the 850 mb level. The ratio, expressed as a 
percentage, is shown in the last column, and is an estimate (probably on the high 
side) of the Sudbury source contribution to the dry deposition at Dorset. 
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Results indicate that the Sudbury sources contributed about 88% of the 
SO2 concentration and 81% of the SO^= concentration of the stratified samples 
during the second operating period. An upper estimate of the Sudbury source 
contribution to the dry SO2 deposition at Dorset site is about 31% for the period 
from July 1981 to March 1982. The Sudbury source contribution to the dry SO^= 
deposition is about 18% for the second period of interest. 

7.2 Railton Air Sample Analysis 

7.2.1 Data Stratification 

The synoptic patterns for the most probable Sudbury smelter-affected 
cases at Railton are similar to those for the Dorset data. These trajectory 
patterns are given in Figure 7-2. 

7.2.2 Data Analysis 

7.2.2.1 Significance of the Sudbury Smelter Contribution 

Since the sample size of the first period (80/81) is too small for the 
detailed data stratification, the estimation of Sudbury smelter contribution for this 
period is not given here. 

Results of the median 2-sample test (see Table 7-1) indicate that the 
SO2 concentrations of the shutdown period (82/83) are significantly lower than 
those of the second operating period (81/82). However, in the cases of SO^ = and N- 
NO3", they are not significantly different. In other words, the Sudbury smelters 
contribute a significant amount of SO2 to Railton, but their contribution of SO^= is 
relatively small. 

7.2.2.2 Estimation of the Sudbury Smelter Contribution to the Dry 
Deposition at Railton 

The mean SO2 concentration and the percentage difference between 

the shutdown and operating periods are given in Table 7-3. The Sudbury smelter 
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contribution to the air SO2 concentration at Railton is about 86%. The Sudbury 
smelter contribution to the dry deposition is estimated by using the approach 
similar to that in section 7.1.2. The estimated Sudbury smelter contribution to the 
SO2 dry deposition at Railton is about 27%. Again, this estimate is thought to be 
on the high side, for reasons already discussed in connection with the Dorset 
analysis. 

7.3. Longwoods Air Sample Analysis 

7.3.1 Data Stratification 

The synoptic patterns for the most possible Sudbury smelter-affected 
cases at Longwoods are given in Figure 6-5. The criteria of selection of the most 
possible Sudbury smelter affected cases are: 

(1) north or northwest of a low type synoptic pattern; 

(2) the 48-hour back-trajectory should at all time be within the north 
sector (45°). 

The criterion used for selecting the possible Sudbury smelter affected 
cases is that the 48-hour back trajectories should at least one time during the 
sampling period be in the north sector (45° sector). 

7.3.2 Data Analysis 

7.3.2.1 Significance of the Sudbury Smelter Contribution 

Since the sample size of the first period (80/81) is not large enough, the 
analysis is not given here. 

Results of the median 2-sample test (see Table 7-1) indicate that only 
SO2 concentrations but not SO^= and NO3- are significantly different between the 
Sudbury smelter operating (81/82) and shutdown (82/83) periods. 
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7.3.2.2 Estimation of the Sudbury Smelter Contribution to the Dry 

Deposition at Longwoods 

The mean air SO2 concentration and the percentage difference between 
the shutdown and operating periods are given in Table 7-3. The Sudbury smelter 
contribution to the air SO2 concentration at Longwoods is about 74% (for the 
stratified data set). The estimated Sudbury smelter contribution to the SO2 dry 
deposition at this site is about 5% (see Table 7-3), a value probably on the high side 
as discussed above. 

7.4 Fernberg Air Sample Analysis 

7.4.1 Data Stratification 

Due to the fact that the sample size is too small after being stratified 
by both the trajectory pattern and the synoptic pattern, only the synoptic pattern 
is used to select the most possible Sudbury smelter affected cases. The criteria of 
selection of the most possible Sudbury affected cases are: that Fernberg should be 
to the northwest of a high pressure system and that the envelope of the 48-hour 
back-trajectories calculated every 6 hours during the 24-hour sampling period 
should include Sudbury. 

7.4.2 Data Analysis 

7.4.2.1 Significance of the Sudbury Smelter Contribution 

Due to the limited sample size of the first operating period (80/81), the 
analysis for this period is not given here. 

The results of the median 2-sample test of the significance of the 
difference between the samples collected during the operating (81/82) and the 
shutdown (82/83) periods are given in Table 7-1. These results indicate that only 
the air SO2 concentrations are significantly different between the two periods. If 
we assume that this difference is mainly due to the Sudbury smelter shutdown (we 
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do not have enough meteorological information to justify this), we may then esti- 
mate the Sudbury smelters contribution to air concentration for the stratified data, 

and to dry deposition at Fernberg. This estimation is given in the next section. 

7.4.2.2. Estimation of the Sudbury Smelter Contribution to Dry Deposition at 

Fernberg 

The result of the calculations is given in Table 7-3. This result 
indicates that Sudbury smelter contributed to the air SO2 concentration at Fern- 
berg is about 57%. By using the approach similar to that for the other sites we 
may estimate the possible Sudbury smelter contribution to the total dry deposition. 
Since samples are available only after October, 1981, an estimate is made for the 
period from October 1, 1981 to March 31, 1982 (see Table 7-2). These results 
indicate that the Sudbury smelters contribute about 4% of the total SO2 dry 
deposition (from October 1, 1981 to March 30, 1982) at Fernberg. 

7.5 Kapuskasing Air Sample Analysis 

7.5.1 Data Stratification 

The synoptic patterns for the possible Sudbury smelter-affected cases 
at Kapuskasing (Ontario Hydro network site) are west or northwest of a high 
pressure system, post-ridge/pre-warm front; warm sector (with southerly flows); 
and pre-north-south oriented cold front (see Figure 7-6). Of these, the case where 
Kapuskasing is west or northwest of a high pressure system is the one which occurs 
most frequently. 

7.5.2 Data Analysis 

7.5.2.1 Significance of the Sudbury Smelter Contribution 

Since only the second and third periods have enough data for analysis, 
the following analysis is for those two periods only. 
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The NO3- concentration of the stratified samples for both the operating 
and shutdown periods are below the analytical detection limit, and no meaningful 
comparison can be made. Results of the median 2-sample test (see Table 7-1) 
indicate that the SO2 and SO^= concentrations during the shutdown period are 
significantly different from those of the 81/82 operating period. These results 
indicate that the Sudbury smelter contribution is significant. 

7.5.2.2 Estimation of the Sudbury Smelter Contribution to the Dry Deposi- 
tion at Kapuskasing 

The mean SO2 and SO^= concentration and the percentage difference 
between the shutdown and operating periods are given in Table 7-3. The Sudbury 
smelter contribution to the SO2 and SO$ a concentrations at Kapuskasing is about 
92 and 66% respectively. The estimated Sudbury smelter contribution to dry 
deposition at Kapuskasing is given in Table 7-3, namely, about 47% and 46%, 
respectively for SO2 and SO^=, which values are again expected to be near upper 
limits, due to the nature of the present analysis. 

7.6 Chalk River Daily Air Sample Analysis 

7.6.1 Data Stratification 

The synoptic pattern and trajectory patterns of the most possible 
Sudbury smelters affected cases given in Figure 6-7. In these cases, the pollutants 
emitted from the Sudbury smelters are transported to the Chalk River 
(Environment Canada APN site) area. 

7.6.2 Criteria of Selection 

The criteria of selection of the most possible Sudbury smelters affected 
cases are: 
1. north or northwest of a high pressure system; 
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2. the 48-hour back trajectories should at all times be within the west sector 
(450). 

To calculate an upper estimate of the contribution of Sudbury smelters 
to dry deposition, cases were used where the 48-hour surface geostrophic back 
trajectories during the 24-hour sampling period were at least one time within the 
west sector between northwest and southwest. Since in general, a westerly 850 
mb-level back trajectory may be associated with a northwesterly to westerly 
surface geostrophic wind trajectory when a cold front is passing through and with a 
southwesterly to westerly surface geostrophic trajectory when a warm front is 
moving away, this criterion should be sufficient for considering the possibility of 
Sudbury smelter contribution from surface to 850 mb-level. 

7.6.3 Data Analysis 

7.6.3.1 Significance of the Sudbury Smelter Contribution 

In general, the air concentrations of the shutdown period are lower than 
those of the operating period. Results of the median 2-sample test (Table 7-1) 
indicate that the SO2 and SO^ = concentrations of the shutdown period are 
significantly lower than those of the operating periods. However, the NO3- 
concentration of the shutdown period is not significantly different from that of the 
operating period. These results indicate that the Sudbury smelter shutdown is the 
major factor responsible for concentration reduction during the shutdown period. 

7.6.3.2 Estimation of the Sudbury Smelter Contribution to the Dry Deposition 
at Chalk River 

The mean SO2 and SCty = concentrations and the percentage difference 
between the shutdown and operating period are given in Table 7-3. The Sudbury 
smelter contributions to the SO2 and SO^= concentration (for the most probable 
Sudbury smelter-affected cases) at Chalk River for the first operating period are 
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about 87% and 55%, and for the second period are about 85% and U7% respectively. 
The estimation of the Sudbury smelter contribution to the dry deposition at Chalk 

River is given in Table 7-3. Results indicate that the Sudbury smelter contribu- 
tions to the dry SO2 and SO^ = deposition for the first operating period are about 
25% and 17% and for the second period are about 26% and 20% respectively. 
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STATISTICAL ANALYSIS SVSTEH 



OHS 



— ►# of data points used 
Period.t*i*- 



8008 
-8844. 



8|03 
8105 
8108 
8111 

8?oa 

-8^0* 



8?08 

8?1L 
8302 



0.0704280 




TABIZ 4-1 SUfARY OF DATA USED IN ANALYSIS 
Hf+ H^ S0 4 3 N-*©3~ Cu ~ 

M£AN1 -MEAN2 *€4*3 M£AN4 m£AN5 



0.0667*52 
0.0**786 7« 
0*0615024 



0.094068 

MH»» 

0.082553 
0.082538 

6.187562 
6.079635 
0.089163 
0.166611 
0.06*981 
6.866387 
6. 061566 



1.327*7 

3^08*27 
2.67350 
3.91586 
1.86526 
1.66727 
1.66737 
4.16*66 
2.69*55 
1.69162 
2.1866* 



0.151679 
-6*63385 1 
6.516367 
6.477636 
6.366928 
6.366666 
6.552663 

6.379669 
6.352766 
6.512666 



0.0011215 
-0.002756*- 
0.0036*9* 
0.00*0550 
0.0021669 
0.0042116 
0.0119957 
6.6017551 
6.6610914 
0.002100* 
0.0017666 



i of missing data points 



HIS*! * HISS2 - MISS3. - MISS4 MlSSS 



1 




."I7TT 



♦Parameters: 1-Bf* f 2=H t + , 3=S0 4 -, 4=r7-!T0 3 -and 5= Cu . (in mg/1 unit) 



Table 4-2 



Statistical Significance at 95X Confidence Level Using 
Median 2-S«*le Test 







July - Sept 




Oct. - Oec 






Jan. - Mar 






1980 vs. 81 


1980 vs 82 1981 vs 82 


1980 vs 81 


1980 vs 82 


1981 vs 82 


1981 vs 8 2 


198 lvs 8 3 


1982 vs 8 3 


«r* 


No 
(.054) 


No Yes 
(.544) (.041) 


No 
(.160) 


No 
(.321) 


No 
(.590) 


No 
(.493) 


No 
(.160) 


No 
(.590) 


< 


No 
(.092) 


No Ves 
(.460) (.001) 


No 
(.388) 


No 
(.422) 


Yes 
(.015) 


No 
(.414) 


Ves 
(.010) 


No 
(1.0) 


»«' 


No 
(.156) 


No Yes 

(.756) (.016) 


No 
(.072) 


NO 
(.214) 


No 
(.590) 


No 
(.072) 


No 
(.327) 


No 
(.281) 


N-N0 3 " 


No 
(.852) 


No No 
(.760) (.218) 


No 
(.695) 


No 
(1.00) 


No 
(.267) 


No 
(.898) 


No 
(.119) 


No 
(.894) 


Cu 


No 
(.104) 


No Yes 
(1.00) (.022) 


No 
(.855) 


Ves 
(.001) 


Yes 

(.000) 


Yes 
(.007) 


Yes 
(.001) 


No 
(.054) 



I 

u> 

U) 

I 



Values In brackets represent probability that the comparison periods' results are not different. 
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TABLE 5.1 



APIOS Cumulative Network Station Identification 



STATIONS 


i 


COLCHESTER 


2 


MERLIN 


3 


PORT STANLEY 


4 


W1LKESP0RT 


5 


ALV1NSTQN 


6 


HUROM PARK 


T 


WATERLOO 


8 


PALMERS TOM 


9 


SHALLOW LAKE 


10 


MILTON 


n 


UXBR1DGE 


12 


COLDWATER 


13 


CAMPBELLFORD 


M 


KALADAR 


15 


SMITH'S FALLS 


16 DALH0US1E MILLS 


17 


GOLDEN LAKE 


18 


WILBERFORCE 


19 


WHITNET 


20 


DORSET 


21 


MCK&LAR 


22 


MATTAWA 



23 


KILLARNEY 


2-4 


BEAR ISLAND 


25 


GOWGANDA 


26 


RAMSEY 


27 


MOONBEAM 


28 


ATTAWAPISKAT 


29 


WINISK 


30 


NAK I NA 


31 


DOR I ON 


32 


QUE T ICO CENTRE 


33 


LAC LA CROIX 


34 


E-L.A. 


35 


EAR FALLS 


36 


PICKLE LAKE 


zr 


BURWASH 


38 


LIVELY 


39 


HANMER 
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Table 6-1 



Results of the Median 2-Sample Test of the Significance of the Difference (precipitation 

concentration) Between the Most Possible Sudbury Affected Samples Collected During 

the Sudbury Shutdown and Operating Periods* 



Site 


Comparison 
Periods 


Parameter 


Test Value 

Z 


Probability 


Significant 




80/81 vs 82/83 


SO^= 
N-NO3- 


-2.13 
-0.98 


0.03 
0.33 


Yes 

No 


Dorset 


81/82 vs 82/83 


N-NO3- 


1.82 
1.01 


0.07 
0.31 


Yes 

No 




80/81 vs 81/82 


so<r 

N-NO3- 


-0.64 
0.64 


0.52 
0.52 


No 
No 




80/81 vs 82/83 


so<r 

N-NO3- 


-1.78 
-0.44 


0.08 
0.65 


Yes 

No 


Railton 


81/82 vs 82/83 


SO/r 
N-NO3- 


-0.44 
-1.78 


0.65 
0.08 


No 
Yes 




80/81 vs 81/82 


soer 
N-NO3- 


1.11 
1.11 


0.27 
0.27 


No 
No 


Longwoods 


80/81 vs 81/82 


so*r 
N-NO3- 


0.89 
0.89 


0.37 
0.37 


No 
No 




80/81 vs 82/83 


SO(t= 

N-NO3- 


1.49 
0.75 


0.14 
0.46 


No 
No 


Chalk River 


81/82 vs 82/83 


SO<,= 

N-N03- 


-2.05 
-0.68 


0.04 
0.49 


Yes 

No 




80/81 vs 81/82 


so^= 
N-N03- 


-0.68 
-0.68 


0.49 
0.49 


No 
No 



* Shutdown period: 3uly 1982 - March 1983. Operating Period: 3uly 1980 - March 1981 
and July 1981 - March 1982. 
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Table 6.2 



Precipitation-depth- weighted Mean and Standard Deviation 

of SO^= and N-NO3- Concentration (mg/1) of the 

Most Possible Sudbury Smelter Affected Cases 



Site 


Period* 




Precipitation 


SO^ : 




N-NO3- 






N 


Total Depth (mm) 


Mean 


S.D. 


Mean 


S.D. 




80-81 


9 


20.1 


1.36 


0.64 


0.24 


0.15 


Dorset 


81 - 82 


3 


11.6 


1.46 


0.41 


0.18 


0.04 




82- 83 


4 


11.4 


0.45 


0.34 


0.14 


0.16 




80-81 


6 


9.7 


1.92 


0.93 


0.84 


0.35 


Railton 


81 -82 


6 


16.8 


2.31 


0.86 


0.89 


0.38 




82- 83 


3 


19.6 


1.69 


0.84 


0.40 


0.19 




80-81 


3 


1.35 


2.50 


0.38 


0.42 


0.13 


Longwoods 


81 -82 


6 


39.7 


2.49 


0.28 


0.27 


0.10 




82-83 


1 


9.9 


1.60 


- 


0.20 


- 




80- 81 


3 


13.2 


6.60 


3.70 


3.03 


2.01 


Chalk River 


81 -82 


5 


49.5 


3.36 


1.49 


1.67 


0.82 




82-83 


3 


17.7 


1.73 


0.43 


1.13 


0.57 



* All periods refer to July of the first year to March of the second year. 



Table 6.3 



Estimated Sudbury Smelter Contribution to Precipitation Concentration and Wet Deposition 



(I) 



OM1M2) 100x(3)/U) 



<4> 



(5)=<3)x<4> (6) ( 7 ) = 100x(5)/(6) 



Site 



Period* 



MPSAC* Average SO*= concen- Excess 
tration(mg I" 1 ) SO*= 

(mgl-l) 



Sudbury con- 
tribution to 

SO%= concen- 
tration (%) 



Possible Sudbury Sudbury con- TotaJ period Sudbury smelters- 
affected cases tribution SO*= Wet SO*= contribution 
Wet deposition deposition S°4 = wet 
(mg m-2) (mg m-2) deposition (%) 



N Operating N Shutdown 



Total precipi- 
tation (mm) 



u> 
-J 

I 



Dorset 80/81 vs 82/83 9 1.36 * 0.45 0.91 

81/82 vs 82/83 3 1.46 * 0.45 1.01 



67 
69 



52 
52 



291 
250 



264.8 
252.5 



2202.7 
2229.3 



12 
II 



Railton 


80/81 vs 82/83 


6 


1.92 


3 


1.69 


0.23 


12 


38 


187.4 


43.1 


2050.0 


2 


Chalk River 


81/82 vs 82/83 


5 


3.36 


3 


1.73 


1.63 


49 


49 


79.8 


227.0 


1955.0 


12 



• All periods refer to July of the first year to March of the second year as indicated. 

♦ MPSAC = Most Possible Sudbury Smelters Affected Cases. 
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Table6-# 



Results of Median 2-Sample Test of the 

Significance of the Difference (Deposition) at Dorset 

Between the Sudbury Smelter Shutdown and Operating Periods* 



Test Test Test Value Probability Significant 

Period Parameter Z >1Z1 



(80-81 and 81-82) vs (82-83) SO 
(80-81) vs (81-82) 
(80-81 and 81-82) vs (82-83) 
(80-81) vs (81-82) 

* All periods refer to July of the first year to March of the second year as indicated. 



SO„ 


-0.365 


0.176 


No 


so 4 


-0.638 


0.523 


No 


N-N0 3 


-0.365 


0.715 


No 


N-N0 3 


0.638 


0.523 


No 
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Table 7.1 

Results of the Median 2-Sample Test of the Significance of the 
Difference of the Air Concentration Between the Shutdown and Operating periods 



Site 


Period* 


Parameter 


Test Value 

Z 


Probability 


Significant 


Dorset 


80/81 vs 82/83 


SO(,= 


-2.77 


0.01 


Yes 






N-NO3- 


AM 


0.07 


Yes 




81/82 vs 82/83 


SCty= 


-2.77 


0.01 


Yes 






N-NO3- 


1.16 


0.24 


No 






S0 2 


-1.66 


0.10 


Yes 




80/81 vs 81/82 


so<r 





1.0 


No 






N-NO3- 


-1.55 


0.12 


No 


Rail ton 


81/82 vs 82/83 


so<r 


-1.08 


0.28 


No 






N-NO3- 


-1.24 


0.22 


No 






so 2 


-2.54 


0.01 


Yes 


Longwoods 


81/82 vs 82/83 


SO$ a 


-Ml 


0.16 


No 






N-N03- 


-Ml 


0.16 


No 






S02 


-3.15 


0.002 


Yes 


Fernberg 


81/82 vs 82/83 


SO^= 


-0.40 


0.68 


No 






N-N03- 


1.01 


0.31 


No 






S02 


2.42 


0.02 


Yes 


Kapuskasing 


81/82 vs 82/83 


SO*" 


2.34 


0.025 


Yes 






so 2 


2.34 


0.025 


Yes 


Chalk River 


80/81 vs 82/83 


so* s 


-2.45 


0.01 


Yes 






N-N03- 


-0.98 


0.33 


No 






so 2 


-2.45 


0.01 


Yes 




81/82 vs 82/83 


so^= 


-2.25 


0.02 


Yes 






N-N03- 


-0.98 


0.33 


No 






S02 


-2.25 


0.02 


Yes 




80/81 vs 81/82 


SCty= 


-0.17 


0.14 


No 






N-N03- 


-0.60 


0.55 


No 






S02 


0.84 


0.40 


No 



* All periods refer to 3uly of the first year and March of the second year as indicated. 
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Table 7.2 

Mean and Standard Deviation of the Air Concentration 
(ug/m 3 ) of the Most Possible Sudbury Smelter Affected Cases 



Site 


Period* 


N 


SO(»= 
Mean 


S.D. 


N 


NO3- 
Mean 


S.D. 


S0 2 

N Mean 


S.D. 


Dorset 


80-81 
81 -82 
82-83 


10 
10 
8 


1.67 
1.73 
0.33 


1.0 
1.1 
0.28 


10 

9 

8 


0.012 
0.016 
0.013 


0.006 
0.018 
0.006 


1 11.23 
9 4.53 
8 0.53 


4.07 
0.51 


Railton 


81 -82 

82-83 


16 
9 


1.71 
1.09 


1.47 
0.96 


16 
10 


0.08 
0.14 


0.10 
0.19 


16 4.66 
9 15.3 


5.74 
13.5 


Longwoods 


81 -82 

82 -83 


11 
7 


2.31 
1.08 


1.46 
0.56 


11 
7 


0.20 
0.10 


0.15 
0.08 


10 3.47 
7 0.90 


3.35 

0.51 


Fernberg 


81 -82 

82- 83 


3 
1 


0.7* 
1.26 


0.39 
0.32 


3 
ft 


0.007 
0.003 


0.006 
0.005 


3 2.65 

4 1.15 


1.13 
0.19 


Kapuskasing 


81 -82 

82- 83 


3 
3 


6.15 
2.08 


3.60 
0.55 


3 

3 


0.80 
0.1 


0.57 

0.06 


J 8.17 
3 0.67 


5.40 
0.12 


Chalk River 


80- 81 
81 - 82 
82-83 


6 

5 

ft 


1.62 
1.38 

0.73 


0.44 
0.81 
0.10 


5 

5 
ft 


0.61 

0.50 
0.27 


0.35 
0.36 
0.13 


5 6.15 
5 5.36 
* 0.83 


5.12 
2.89 
0.59 



• All periods refer to 3uly of the first year and March of the second year as indicated. 











Estimatec 


1 Sudbury Smelter 


Table 7.3 
Contribution to Dry Deposii 


tion and Air Concentration 














(1) 




(2) 


(3M1M2) 


I00x(3)/(l) 


(4) 


(5)=(3)x(4) 


(6) 


(7)= 


:100x(5)/(6) 


Site 


Period* 


Parameter 


MPSAC* Average 
Operating 


Concentration 
Shutdown 


Excess 
(ug/m*) 


Sudbury contri- 
bution to air 
concentrator 


Possible 
Sudbury affec- 

t ted cases 


Estimated 
Sudbury con- 
trib.(ug/m3) 


Total concen- 
tration of the 

period (ug/m3) 


Sudbury Smelter 
contribution (%) 


Dorset 


80-81 
81 -82 


so 2 
so«= 
so 2 

SO* 3 


N 

(4) 
(10) 
(9) 
(10) 


1.67 
4.53 
1.73 


N 

(10) 
(8) 
(10) 
(8) 


0.33 
0.53 
0.33 


4.0* • 
1.34 
4.0 
1.40 


80 
88 
81 


N 

170 
170 
127 
127 


680 
228 
508 
177 


1842 
810 
1616 
969 




37 
28 
31 
18 


Rail ton 


81 -82 


so 2 


(16) 


4.66 


(8) 


0.66 


4.0 


86 


127 


508 


1902 




27 


Longwoods 


81 -82 


so 2 


(10) 


3.47 


(7) 


0.90 


2.57 


74 


68 


175 


3250 




5 


Fernberg 


81 -82 


so 2 


(3) 


2.65 


(4) 


1.15 


1.5 


57 


13 


19.5 


492 




4 


Kapuskasing 


81 -82 


so 2 
so*= 


(3) 
(3) 


8.17 
6.15 


(3) 

(3) 


0.67 
2.08 


7.5 

4.07 


92 
66 


46 
46 


345 
187 


737 
406 




47 
46 


Chalk River 


80-81 
81-82 


so 2 
so»= 
so 2 
so*= 


(6) 
(6) 
(4) 
(4) 


6.15 
1.62 
5.36 
1.38 


(4) 
(4) 
(*) 
(4) 


0.83 
0.73 
0.83 
0.73 


5.32 
0.89 
4.53 
0.65 


87 
55 

85 

47 


121 
121 
142 
142 


644 
108 
643 
92.3 


2572 
625 

2468 
468 




25 

17 
26 

20 



I 

H 

1 



• MPSAC = most possible Sudbury affected case 

• All periods refer to July of the first year to March of the second year as indicated. 

• • The 80/81 excess concentration is assumed to be the same as that of 81/82. 
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FIGURE 2 - 1 : STATION LOCATION MAP 

CUMULATIVE PRECIPITATION 
MONITORING NETVCRK 



« 



35 



6 



30 



27 



31 



32 



*• 



«*£> 



ii 



ONTARIO 


t t> m m m 


km 


^•■^ Enmro#vn«nf 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



COLCHESTER 

MERLIN 

PORT STANLEY 

WILKESPORT 

ALVINSTON 

HURON PARK 

WATERLOO 

PALMERSTON 

SHALLOW LAKE 

MILTON 



11 - UXBRIDGE 21 

12 - COLDUATER 22 

13 - CAMPBELLFORO 23 

14 - KALAOAR 24 

15 - SMITH'S FALLS 25 

16 - OALHOUSIE MILLS 26 

17 - GOLDEN LAKE 27 

18 - WILBERFORCE 28 

19 - WHITNEY 29 

20 - DORSET 30 



MCKELLAR 

MATTAMA 

KILLARNEY 

BEAR ISLAND 

GOUGANOA 

RAMSEY 

MOONBEAM 

ATTAWAP ISKAT 

WINISK 

NAKINA 



31 - DOR I ON 

32 - OUETICO CENTRE 

33 - LAC LA CROIX 

34 - EXP. LAKES AREA 

35 - EAR FALLS 

36 - PICKLE LAKE 
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FIGURE 5-1 IS04 CONCENTRATION (MG/l) FROM JULY 1980 TO MARCH 1981 
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FIGURE 5-2:S04 CONCENTRATION (MG/U FROM JULY 1981 TO MARCH 1982 
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FIGURE 5-3:S04 CONCENTRATION (MG/L) FROM JULY 1982 TO MARCH 1983 
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F1GURE5-4*AT10 OF S04 CONCENTRATION BETWEEN 82-83 AND 80-81 
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FIGURE 5-5,RAT 10 OF S04 CONCENTRATION BETWEEN 82-83 AND 81-82 
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FIGURE 5-6:N-N03 CONCENTRATION (MG»100/U FROM JULY 1980 TO MAR 1981 



'6' 



-55- 



u» #4* t2 § to* ■•• !«• »4* »a» io' t%* re* rv tv 




S i ■ so 



n* to* ••■ •«• i4« ia« io* r»« 

FIGURE 5-7: M-M03 CONCENTRATION (MG«100/U FROM JULY 1981 TO MAR 1982 
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FIGURE 5-8 ; N-N03 CONCENTRATION (MG»100/U FROM JULY 1982 TO MAR 1983 
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FIGURE5-9,RAT10 OF N-N03 CONCENTRATION BETWEEN 82-83 AND 80-81 
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FIGURE 5-10: RATIO OF N-N03 CONCENTRATION BETWEEN 82-83 AND 81-82 
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FIGURE 5-11: S04/N-N03 RATIO FROM JULY 1980 TO MAR 1981 
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FIGURE 5-12: S04/N-N03 RATIO FROM JULY 198l TO MAR 1 982 
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FIGURE 5-13: S04/N-N03 RATIO FROM JULY 1982 TO MAR 1983 
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FIGURE 5-14: RATIO OF S04/N-N03 RATIO* OF 82-83 AND 80-81 
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FIGURE 5-15: RATIO OF SO4/N-NO3 RATIOS OF 82-83 AND 81-82 
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FIGURE 5-16 : CU CONCENTRATION (UG/L) FROM JULY 1980 TO MARCH 1981 



-65- 



rr TV TV T2* 




•a* to* »•• •«• «4« »a* to* 

FIGURE 5-17 : CU CONCENTRATION (UG/L) FROM JULY 1981 TO MARCH 1982 
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FIGURE 5-18 I CU CONCENTRATION (UG/U FROM JULY 1982 TO MARCH 1983 
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FIGURE 5-19 : S04 DEPOSITION (MG«10/M««2) FROM JULY 1980 TO MAR 1981 



7%' 



-68- 



Mi 94' 92* 80* 88* 86* 84* 82' 80* 78' TV 74' 




6-90 



92* 80* 88* 86* 84* 82* 80* 

FIGURE 5-20 : S04 DEPOSITION (MG«tO/M««2) FROM JULY 1981 TO MAR 1982 



-69- 



H* «4» «* 90" !•• »«• »4» »2* 80 % /•' TV 



7 4' 




2.10 



da* «o* ••• M* |4» «• so* n« 

FIGURE 5-21 : S04 DEPOSITION <MG«10/M««2) FROM JULY 1982 TO MAR 1983 
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FIGURE 5-22:- RATIO OF S04 DEPOSITION BETWEEN 82-83 AND 80-81 
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FIGURE 5-23 : RATIO OF S04 DEPOSITION BETWEEN 82-83 AND 81-82 
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FIGURE 5-24: N-N03 DEPOSITION tHG/tf**2J FROM JULY 1980 tq MAR 1981 
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FIGURE 5-25: N-N03 DEPOSITION (MG/M««2) FROM JULY 1 98 1 TO MAR 1982 
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FIGURE 5-26: N-N03 DEPOSITION (MG/M««2) FROM JULY 1982 TO MAR 1983 
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FIGURE 5-27: RATIO OF N-N03 DEPOSITION BETWEEN 82-83 AND 80-81 
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la) Post-cold Front Pre-high Pressure Ridge 




lb) Cold-front Moving Through 




lc) Southwest of a Low Pressure System 



figure 6-1. Synoptic Patterns of Possible Sudbury Smelters 
Contributed Cases. 
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Figure 6-? 



Precipitation Concentration of the Most 
Possible Sudbury Smelters Affected Cases 
at Dorset. ( ® denotes ? points.) 
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Figure 6-2. (Cont'd) 
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Figure 6-3. Deposition of the Most Possible Sudbury 
Smelters Affected Cases at Dorset 
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Figure 6-3. (Cont'd) 
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Figure 6-4. Synoptic Pattern of the Most Possible 
Sudbury Smelters Affected Cases at 
Railton. 
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Figure 6-5. Synoptic Pattern of the Most Possible Sudbury 
Smelters Affected Cases at Longwoods. 
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a. Northwest of a high pressure system 




b. Southwest of a high pressure system 




c. Northeast of a Low oressure system 



Figure 6-6. Synoptic Patterns of the Most Possible 

Sudbury Smelters Affected Cases a.t Fernberg. 
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Figure 6-7. 



Synoptic Pattern - of the Most Possible Sudbury 
Smelters Affected Cases at Chalk River. 
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Figure 7-1. 



^6-hour Back-trajectory Patterns of the Most 
Possible Sudbury Smelters Affected Cases at. 
Dorset . 
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Figure 7-2 



48-hour Back -trajectory Patterns of the Second Most 
Possible Sudbury Smelters Affected Caner, 
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Figure 7-3. Ground-level Air SO Concentration of the Most 

Possible Sudbury Smelters Affected Cases at Dorset. 
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Figure 7-4. Cround-level Air S0 ? Concentration of the moat 

Possible Sudbury Smelters Affected Cases at Dorset. 
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Figure 7-5. Ground-level Air NO Concentration of the Most 
Possible Sudbury Smelters Affected Cases at 



Dorset. 
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a. West or Northwest of 
a high 



b. Post-ridge/pre-warm front 





c. Warm sector 

(with southerly flows ) 



d. Pre-cold front 

(north-south oriented) 



Figure 7.6. 



Synoptic Patterns of the Most Possible Sudbury 
Smelters Affected Cases at Kapuskasing. 
(Receptor area is hatched. Source area is to 
the S of the receptor.) 
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APPENDIX 1 

Sudbury Shutdown Study Data Stratification Method for Daily 
Precipitation Samples at the Dorset site 

This appendix discusses the data stratification method used for daily precipi- 
tation samples, and the related synoptic conditions. 

1. Background 

Due to the fact that INCO emission sources are not the only important 
sources and that the samples collected are subject to the effects of complicated 
meteorological processes, the data have to be properly stratified to take into 
account the important meteorological processes. 

The present stratification scheme attempts to group together data which 
were obtained under similar (or the same) meteorological (transportation, dilution, 
stagnation, wet/dry deposition) and chemical (transformation, oxidation) condi- 
tions. Traditionally, data stratification is done by stratifying data only according 
to trajectory sectors, to account for the different source strength in different 
compass direction. However, this is not sufficient, since a trajectory sector is not 
associated with a unique synoptic condition. In other words, a simple trajectory 
analysis does not fully account for the meteorological variability. As a conse- 
quence, trajectory sector statistics vary every year even when the emission rates 
of all sources within a given sector are constant. Besides, a given trajectory 
sector, unless very carefully selected, may not relate to a unique combination of 
sources, hence, the stratified data still may have a wide range of variation. 

Data are stratified by another frequently used method according to synoptic 
patterns, since similar patterns are associated with similar meteorological condi- 
tions. However, this method cannot easily take the emission sources into account, 
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since the synoptic pattern changes with time at a specific site. The variation of 
the affecting sources in a synoptic pattern makes the statistics of a synoptic 
pattern vary between years, and the data associated with a specific synoptic 
pattern has a wide range. 

For these reasons, the method used here was developed. However, it is first 
necessary to understand clearly the meteorological situations when precipitation 
at the receptor is occurring. In the following section, the important meteoro- 
logical considerations related to data stratification are given. 

2. Meteorological Considerations. 

Precipitation may be classified as frontal and non-frontal types. For 
the non-frontal type, the use of single level trajectory (surface or 850 mb level) 
may be sufficient to distinguish between sources. However, the frontal type needs 
special considerations. 

For trajectory analysis to be useful, the trajectories must correctly indicate 
the origin of the polluted air masses which have brought about most of the 
deposition at the receptor. In frontal situations, this may not always be the case, 
since the flows that move up along the front into the cloud neither follow a 
constant pressure surface nor a constant height surface (Figure A-l). Surface 
geostrophic trajectories indicate the correct surface air parcel movement only if 
the calculated trajectories do not intercept the surface front (since no air parcel 
moves across the front). In other words, the surface geostrophic trajectories may 
only be used to estimate the air parcel origin when they are in a cool sector (before 
warm front), a warm sector (after warm front, before cold front) and in a cold 
sector (after cold front). 

For example, the Sudbury emissions may contribute to Dorset precipitation 
samples when the samples are collected in a cold sector (this includes three types 
of patterns: (1) post cold front: (2) cold front passage; (3) precipitation to the 
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southwest of a low pressure system). In these cases surface geostrophic trajec- 
tories may be used to stratify the Dorset precipitation samples, since the Sudbury 

sources' emissions are transported mainly by cold sector flows (below cold front) to 
the Dorset region. However, especially in the case of frontal passages, trajectories 
at different vertical levels have to be examined due to wind shear and the elevated 
nature of the Sudbury emissions. If an attempt is made to get around this problem 
by widening the sector for which Sudbury has a potential influence (and using only 
geostrophic trajectories, for example), then the variability due to meteorological 
factors is increased. In this present analysis, we looked at the surface geostrophic 
trajectories and the synoptic weather maps to decrease meteorological variability. 

Synoptic patterns have been classified into the following types defined by 
Heidorn (Water, Air and Soil Pollution, U_, 225-235 (1979)). 

(1) post cold front/pre-high pressure ridge 

(2) high pressure ridge 

(3) post-ridge/pre-warm front 

(4) warm front 

(5) cyclonic 

(6) cold front 

(7) stationary front 

(8) weak pressure gradient 

This classification which was originally developed for a climate study, was 
here extended for the present study, taking into account extra factors such as 
stagnation, transporting flow direction and rate, and precipitation location with 
respect to the synoptic pattern. 

This resulted in the following classification scheme (Figure A-2). 
1. Post-cold fron/pre-high pressure ridge 

(a) northwest to the low pressure system 

(b) southwest to the low pressure system 
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2. Warm frontal and warm sector 

This type is ususlly associated with a stagnation high pressure system (located 
in southeast U.S. continent) which moves slowly eastward into the ocean. The 
southerly flow brings heavily polluted air into the warm frontal clouds. This is a 
typical episode type synoptic patttern. 

3. Post-ridge/pre-warm front 
k. Warm front 

5. Warm sector 

6. Cold front 

7. East-west oriented cold front or quasi-stationary front 

(a) north of the front 

(b) south of the front 

8. Northwest of a high pressure system 

9. Southeast of a high presure system 

10. East of a high/ warm frontal 

11. Cyclone 

(a) NE of cyclone 

(b) NW of cyclone 

(c) SE of cyclone 

(d) SW of cyclone 

(e) centre of cyclone 

12. Occluded front 

(a) pre-occluded front 

(b) post-occuluded front 

13. Others. 

For the precipitation Sudbury impact study, category 1 above was examined 
in detail with the Dorset data, for example, because it is the situation where there 
is the highest probability that a Sudbury smelter effect will be seen at Dorset. 
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Since the synoptic system is not stationary, it is necessary to classify the 
sample by all the possible combinations of the above classes. In the shutdown 
study, samples corresponding to mixed synoptic patterns should not be used in 
sample comparisons. 

From the above considerations, it is clear that we need both trajectory and 
synoptic pattern to stratify data instead of only stratifying according to trajectory 
or synoptic pattern alone . Using both trajectory sector and synoptic pattern to 
stratify the precipitation data, between-year variation of the average sample 
concentration can be minimized. For example, the average concentration of the 
class of post-cold f ront/pre-high ridge class and with trajectories in the NW sector 
did not change significantly between years if the strength of the emission sources 
remain constant. 

It has to be mentioned here that the use of trajectory sector and synoptic 
pattern to stratify data does not fully account for meteorological variability, since 
the above synoptic pattern classification does not fully account for the storm track 
or geometric shape of the synoptic pattern. The use of a trajectory sector cannot 
overcome this problem, since trajectories can have different patterns within a 
sector. In other words, the data after being stratified by trajectory sector and 
synoptic pattern, still may show the within group variation; however, the between 
group mean variation is minimized. If the sample size is large enough, we may 
stratify data by using trajectory pattern and synoptic pattern instead of using 
trajectory sector. This can minimize the within group data range and provide more 
specific statistics. 
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a. Flow on horizontal plane 
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b. Vertical cross section. 




c. 3 dimensional 



Figure A-l :Flow structure at frontal zone 
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Class 1. Post-cold front /pre-high pressure ridge. 



Figure A-2: Synoptic pattern classes. 

(Hatched area denotes receptor region. ) 
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Class 2. Warm frontal and warm sector. 




Class 3. post-ridge/pre-warm front. 
Figure A-2: (Cont'd) . 
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Class 



4. Warm front 




Class 5. Warm sector 
Figure A-2: (Cont'd) 
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Class 6. Cold front 




Class 7. East-west oriented cold front or quasi-stationary 
front . 



Figure A-2:(Cont'd) 
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Class 8. Northwest of a high pressure system 




Class 9. Southeast of a high pressure system 
Figure A-2 (Cont'd) 
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Claas 10. East of 
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Class 11. Cyclone 
Figure A-2 (Cont'd) 



-103- 




Clasa 12, Occluded front 
Figure A-2 (Cont'd) 
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